Quantum Zeno suppression of three-body losses in Bose-Einstein condensates 
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We study the possibility of suppressing three-body losses in atomic Bose-Einstein condensates via 
the quantum Zeno effect, which means the delay of quantum evolution by frequent measurements. 
It turns out that this requires very fast measurements with the rate being determined by the spatial 
structure of the three-body form factor, i.e., the point interaction approximation <5 3 (r — r') is not 
adequate. Since the molecular binding energy Et provides a natural limit for the measurement rate, 
this suppression mechanism can only work if the form factor possesses certain special properties. 
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To the best of our knowledge, Bose-Einstein conden- 
sates of trapped dilute alkali gases are the coldest systems 
in the Universe Apart from their low temperature, 
atomic Bose-Einstein condensates are an ideal laboratory 
for studying many-body quantum effects @, [3] because 
they can be well manipulated and controlled experimen- 
tally by using lasers and magnetic fields [4]. However, 
since the gaseous state is not the real ground state of 
these atoms, inelastic collisions due to three-body recom- 
bination occur [5(. In such an exothermic process, three 
atoms collide to form a dimer and an unbound atom, 
which carries away the excess energy and momentum. 
The binding energy thus released is converted into the 
kinetic energies of the colliding partners (molecule and 
third atom) which leave the trap. As a result of this 
process, the condensate is susceptible to loss and heat- 
ing 0, 0] , thereby inevitably limiting its lifetime and the 
minimum achievable temperature. Consequently, one is 
lead to the question as to how one could suppress this 
undesired loss mechanism. 

In this work we address this important question us- 
ing the idea of quantum Zeno effect which is basically 
the delay or even inhibition of the time evolution of a 
quantum state by performing frequent measurements ]g- 
|lfjj. As an example, consider a simple two- level system 
undergoing Rabi oscillations of frequency lj. When the 
initial state gets rotated under the action of the total 
Hamiltonian H , the state evolves to 
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Without measurements, the probability to find the state 
in |V>i) after time t is P(t) = \(ipi\e- iHt \ipi}\ 2 = cos 2 wt. 
When N measurements are made during this time (after 
equidistant time intervals), the survival probability reads 
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where 1/r 2 = (V>i|# 2 |i/>i) = 1 1] - For large N, the 
survival probability P/v(i) approaches one, i.e., under 
continuous measurements the time evolution of the sys- 
tem is prevented. Thus, quoting 11| Zeno's quantum- 



mechanical arrow (the wave function) , sped by the Hamil- 
tonian, does not move, if it is continuously observed. 

The above example demonstrates the suppression of 
transitions between two discrete levels in a quantum sys- 
tem. However, the aforementioned three-body losses cor- 
respond to a decay into a continuum of states. Thus, we 
have to study whether and when the quantum Zeno ef- 
fect occurs in such cases. To this end, we consider 12 1 
the following model Hamiltonian (h = 1) 

H = E \g)(g\ + J dk [E k \k)(k\ + (g k \g)(k\ + h.c.)] , (3) 

where (g\g) = 1, (k\k') = S(k - k f ), and (g\k') = 0. Here 
\g) denotes an initial quasi-bound state (i.e., a resonance) 
with energy Eq coupled to a continuum of states |fc) with 
energies E k via the coupling strength . Similar Hamil- 
tonians have been used to study, for example, a two-level 
atom interacting with radiation field [1 31 ) . 

To the lowest order in perturbation g k , the total decay 
probability P d (t) = j dk\(k\e~ iAt \g)\ 2 , takes the form 
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where we have used \g(E)\ 2 = \gu\ 2 (dEk / dk)^ 1 . For 
short times, we obtain a quadratic growth of the decay 
probability Pd(t) cx t 2 similar to the previous example in 
Eq. ([T]). Consequently, the time evolution can be delayed 
via the quantum Zeno effect according to Eq. @. For 
longer times (but still within the region of perturbation 
theory), the decay probability grows linearly Pd(t) oc t 
and we recover Fermi's golden rule. In this regime, mea- 
surements would not slow down the time evolution signif- 
icantly. Therefore, the crossover time T* when the transi- 
tion between these two different regimes takes place [III , 
i.e., P d {t < T«) oc t 2 and P d {t > T») cx t, determines 
how fast one has to measure in order to slow down the 
decay via the quantum Zeno effect. We may estimate the 
crossover time T* via calculating the intersection point 
of the two asymptotic expressions Pd{t -C T») oc t 2 and 
Pd(t > T*) oc t which yields 
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This is reminiscent of Heisenberg's energy-time uncer- 
tainty relation and implies that one has to measure be- 
fore the energy conservation begins to play a role in order 
to observe the quantum Zeno effect. 

Alternatively, using the wave-packet description one 
may write T» in the following way 

l = 0(v S[0Up Ak) = 0{^). (6) 

Here Afc and AE are the widths of \g k \ 2 and \g{E)\ 2 , 
respectively, and w gr0 up = dE^jdk is the group velocity 
(such that Vgj-oupA/c rs AE). For early times, the wave- 



packet created by the decay of the initial quasi-bound 
state \g) into the continuum of states \k) possesses the 
width AA; in momentum space and thus it spreads over 
Ar = 0(1/ Ak) in position space. Thus, T* is the period 
during which this wave-packet of spatial size Ar needs to 
move away with its group velocity i> gr0 up- In other words, 
as long as the wave-packet is close to its initial position 
(and thus has large overlap with the original quasi-bound 
state \g)), its decay can be suppressed via the quantum 
Zeno effect - once the wave-packet has crossed its own 
dimensions and moved away, measuring does not change 
much anymore. 



Studying the simple toy model in Eq. (|3J|, we found that the quantum Zeno effect may also occur in the decay to 
a continuum of states, but only if measurements are performed fast enough - with the characteristic time T* being 
determined by the form factor |<?fc| 2 . Now let us study a more realistic model Hamiltonian for three-body losses 
in Bosc-Einstcin condensates. For simplicity, we restrict our consideration to one species of atoms with mass m 
(described by the field operator xjj) coupled to one species of molecules with mass 2m (and the field operator (f>) 
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Here Eb is the negative molecular binding energy and W(r — r') is the atom- molecule coupling. Since the characteristic 
length scales of molecule formation are typically much smaller than the trap size, we assume a homogeneous condensate 
with density p and omit the trap potential V(r). Within standard perturbation theory (i.e., expansion into powers of 
the coupling strength W) we recover the usual Feynman diagrams depicted in Fig. [TJ Three-body losses are described 
by processes of third (or higher) order (i.e., Fig. [TJ;) and the decay probability reads to lowest order in W 
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Here W(k) = W(k) is the Fourier transform of W(r — r') and we have abbreviated lo\ — Eb + 3Ek/2, u 2 = Eb — Ek/2, 
CJ3 = 3Ek/2, and E k = k 2 /(2m). Initially, at time t = 0, we consider a situation \^f(t = 0)) = \^q) where all the atoms 
have Bose condensed ip |^o) = \[f> |^o) an d there are no molecules <fi \^o) = 0. Of course, this is a strong simplification 
(assuming perfect measurements and weak depletion etc.) but the main results should remain unaffected. 



The expression for the time dependence of the decay 
probability in Eq. (JS) allows us to study the occurrence of 
the quantum Zeno effect. For small t, a Taylor expansion 
of the integrand in Eq. ((5|) in powers of t yields 
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Thus, the decay probability increases as t 6 initially. This 
departure from the conventional t 2 behavior obtained 
earlier is attributed to the three-body loss as a third- 
order process. In this regime, frequent measurements 
slow down the decay significantly. 

For large t, on the other hand, we recover the linear 
time dependence analogous to Fermi's golden rule: In 
this limit, the points of stationary phase in the integrand 
dominate the integral in Eq. ([5]). Since the binding en- 
ergy Eb is negative, we have L02 < 0, W3 > and thus 



> o>2 for all k > 0. Consequently, the only point of 
stationary phase (which is also the only point where one 
of the denominators vanishes) is located at u>i = and 
(as one would expect) corresponds to energy conservation 
3Ek/2 = —Eb- Thus we derive the asymptotic (large t) 
behavior of Eq. © in complete analogy to Eq. ((U) : For 
large values oft, only the term containing sin 2 (u;i£/2)/a; 2 
exhibits a peak at u>i = 0, of height cx t 2 and width 
<xl/t. Thus the significant contribution to the integral 
arises from the neighborhood of this peak. Performing 
the integration we get (for large t) 
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where the value of fco = (^mEb/3) 1 ^ 2 is dictated by en- 
ergy conservation. 

The crossover time T* between the two regimes © 
and (|10p can again be estimated by equating these two 
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asymptotic expressions for Pd{t) 
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The integral in the equation above scales with the cubed 
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width fc?„ of the function VF(k) in fc-space (relative to its 



value at k and 2k Q , respectively). Thus, the crossover 
time T* can be estimated as 
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Similar to the previous example in Eq. (J3J, this expres- 
sion facilitates an intuitive interpretation - even though 
it is more qualitative than the quantitative picture de- 
scribed after Eq. ([6]). At early times, the spatial size Ar 
of the molecular wave-packet scales with Ar = 0(1/ kw) 
while its group velocity is determined by ko, he., the bind- 
ing energy Et,. Hence T* grows for both kw and/or ko 
shrinking. 




FIG. 1: Feynman diagrams for the Hamiltonian in Eq. (0. 
Single lines denote atoms and double lines are molecules while 
the black dot is the interaction vertex. The left diagram 
a) shows the first-order process of direct molecule formation 
which is forbidden by energy conservation. The middle di- 
agram b) corresponds to atom-atom scattering via an inter- 
mediate ("virtual") molecular state, which is a second-order 
effect. Finally, the right diagram c) combines the processes 
a) and b) and shows the actual third-order three-body recom- 
bination process resulting in the formation of a real molecule 
plus an atom (carrying away the excess energy and momen- 
tum). 

In units of the inverse binding energy Ef,, the crossover 
time T* scales with the ratio ko/kw to the power 3/5. In 
order to achieve a significant delay of the time evolution 
via the quantum Zeno effect, we would have to measure, 
i.e., detect the formed molecule, faster than this time 
scale . The presence of the molecule could be detected 
via driving a light induced transition between molecular 
bound states which is off-resonant for all atomic states 
(see also [HI])- Then the energy-time uncertainty relation 
implies that this measurement cannot be done faster than 
the inverse binding energy E\, . Even if we could imagine a 
detection method for which this energy-time uncertainty 



bound does not apply, it would be unwise to measure 
faster than the inverse binding energy Ef,: For such mea- 
surements, the direct two-body recombination process in 
Fig. [1^, which is forbidden by energy conservation, would 
start to contribute. Since the rate of this first-order pro- 
cess (if allowed) is typically much larger than that of 
the third-order effect in Fig. [TJ;, the condensate would 
become even more unstable. Therefore, stabilizing the 
Bose-Einstein condensate via the quantum Zeno effect is 
only possible if T„E b > 1. With Eq. (12} we sec that 
this requires ko ^> kw- So if the structure of W(\s) is 
basically determined by just one length scale (e.g., the s- 
wave scattering length a s ) as in Fig. [2^,, we cannot stabi- 
lize the Bose-Einstein condensate via the quantum Zeno 
effect. However, if VF(k) has a more complicated, e.g., 
spiky, structure as in Fig. \23p, we could have ko l$> kw 
and thus it could be possible. 




FIG. 2: Sketch of two possible forms for the atom-molecule 
coupling W as a function of k illustrating the two character- 
istic cases: In the left picture a), W has no sub-structure and 
thus ko and kw are of the same order of magnitude kw ~ ko, 
i.e., the quantum Zeno effect is not applicable. (Note that 
kw 2> ko would still be possible in this case - but in this 
limit, the quantum Zeno effect does also not apply.) In the 
right picture b), the function W has more structure and thus 
we have kw <fC ko, which enables us to apply the quantum 
Zeno effect. 

In summary, we studied the possibility of suppress- 
ing three-body recombinations in dilute atomic Bose- 
Einstein condensates via the quantum Zeno effect. This 
fascinating phenomenon corresponds to the delay or even 
inhibition of quantum evolution by repeated measure- 
ments if they are fast enough. Addressing the ques- 
tion of how fast is fast enough, we first see that the 
probability for a positive measurement outcome (saying 
that the system has actually evolved) must always be 
small. In addition to this trivial bound, another time- 
scale arises in the decay of a meta-stable state into a 
continuum: In such a set-up, we obtain a quadratic (or 
higher polynomial) growth of the decay probability at 
early times ("quantum Zeno regime") and a linearly in- 
creasing decay probability at later times ( "Fermi golden 
rule regime"). Obviously, measurements can only slow 
down the evolution during the initial period (i.e., the 
"quantum Zeno regime" ) . For the model Hamiltonian in 
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Eq. ([3]), describing the tunnelling of a particle out of a 
potential well, for example, the crossover time T* sepa- 
rating the two regimes [20] is given by the width of the 
corresponding wave-packet in comparison with its group 
velocity. In the case of three-body recombinations in 
dilute atomic Bose-Einstein condensates, this crossover 
time T* is determined by the spatial structure of the 
three-body form factor W(r — r') in Eq. ([7]). As a result, 
the often used [HI, [n| point interaction approximation 
W(r — r') oc 5 3 (r — r') is not adequate for this question. 
Since the molecular binding energy £& provides a natural 
limit for the measurement rate, this suppression mecha- 
nism can only work if the typical width kw of the form 
factor W is much smaller than characteristic momentum 
of the final molecule/atom k = (AmEb/3) 1 / 2 [17 1. 

It might be illustrative to compare our approach with 
the work [13 ] in which the inhibition of three-body re- 
combination via resonant 2ir laser pulses was studied. 
First of all, it should be noted that these 2ir laser pulses 
considered in 14{ do not correspond to an actual mea- 
surement but induce a phase shift. Therefore, the mech- 



anism studied in LJJ is not the quantum Zeno effect, but 
corresponds to passive error correction techniques known 
in quantum information theory under the names "spin- 
echo" or "bang-bang" method, see, e.g., [lit ]. Neverthe- 
less, both phenomena are based on the crucial difference 
between summing amplitudes and probabilities in quan- 
tum theory. Therefore, the repetition rate of the 2n laser 
pulses considered in 14] should lie in the "quantum Zeno 
regime" in order to induce a significant slow-down (in the 
"Fermi golden rule regime" they would have very little 
effect), i.e., similar temporal constraints should apply in 
both cases. 
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